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■ Abstract 

^ ' We construct an effective Hamiltonian via Monte Carlo from a given action. This Hamilto- 
^ , nian describes physics in the low energy regime. We test it by computing spectrum, wave 
d ' functions and thermodynamical observables (average energy and specific heat) for the free 

system and the harmonic oscillator. The method is shown to work also for other local 

potentials. 
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1 Motivation 



The motivation for constructing a Monte Carlo Hamiltonian comes from different direc- 
tions. 

(i) The renormahzation group a la Kadanoff-Wilson [Q] aims to construct a renormalized 
Hamiltonian, which describes physics at a critical point, based on the assumption of scale 
invariance. Such Hamiltonian is supposed to have much less degrees of freedom than the 
original Hamiltonian. A recent further development of those ideas is White's density matrix 
renormahzation group technique 0. Our goal is similar to the above in the sense that we 
aim at an effective Hamiltonian, which has "less" degrees of freedom than the "original" 
Hamiltonian. But it differs in describing physics in the low energy domain instead of doing 
so at the critical point. 

(ii) When one tries to solve field theory in the Hamiltonian formulation, the standard way 
to proceed is by constructing a Fock space, parametrized by some high momentum cut-off 
and some occupation number cut-off (Tamm-Dancoff approximation). When increasing 
those parameters, which means increasing the upper bound of the energy, then typically 
the density of states increases in an exponential manner, which renders the system beyond 
any control. In contrast to that, the Monte Carlo Hamiltonian is governed by a "small" 
number of low-energy degrees of freedom and the spectral density decreases with increasing 
energy. 

(iii) The enormous success of lattice field theory over the last quarter of the century is 
certainly due to the fact that the Monte Carlo method with importance sampling is an 
excellent technique to solve high dimensional (and even "infinite" dimensional) integrals. 
Conventionally, one computes a transition amplitude of an operator and evaluates it nu- 
merically via Monte Carlo (e.g. Metropolis algorithm 0), 

^ I[dx] exp(-^g[x]) 
J[dx] exp{-^S[x]) 

- jrT.o[c]. (1) 

Here C stands for a path configuration drawn from the distribution P[x] = ^ exp(— ^^[x]). 
The vitue of the Monte Carlo method lies in the property of yielding very good numerical 
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results. E.g., solving a field theory model on a lattice of size 20'^ and measuring the 
observable from a number of configurations A^^^ in the order of a few hundred typically 
yields results with statistical errors in the order of a few percent. In this way it has been 
possible to determine low lying baryon and meson masses quite precisely [Q. 

On the other hand, one can express a transition amplitude in imaginary time via the 
Hamiltonian 

oo 

= < ^M^n> e-^"^/" < En\Xin> 

n=l 

N 

= E < ^MK" > e-^^"'^/' < Et^^lx.^ > . (2) 

In the last two lines we have approximated the Hamiltonian H by an effective Hamiltonian 
Heff, which has less degrees of freedom, e.g., it has only eigenstates. The idea of the 
Monte Carlo Hamiltonian is that an effective Hamiltonian can be found via use of Monte 
Carlo, such that transition amplitudes become a finite sum over eigenstates, where N is 
in the order of magnitude of N^, i.e. the number of equilibrium configurations, sufficient 
to closely approximate the path integral of Eq. (0) . 

One might ask: What is the virtue of such a Hamiltonian? A list of physics problems, 
where progress has been slow with conventional methods including standard lattice tech- 
niques, and where such a Hamiltonian might bring progress are the following topics: 

- Non-perturbative computation of cross sections and decay amplitudes in many-body 
systems]^. 

- Low-lying but excited states of the hadronic spectrum and the related question of quan- 
tum chaos in such a system. 

- Hadron wave functions and the related question of hadron structure functions, in particu- 
lar for small xb and Q^. The Hamiltonian formulation is suited to compute wave functions, 
which is quite difficult in the Lagrangian lattice formulation. 

- Finite temperature and in particular finite density in baryonic matter. This is crucial 
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for the quark-gluon plasma phase transition, the physics of neutron stars and cosmology. 
The Hamiltonian formulation is suited to compute the mean value of the energy (average 
energy). This is difficult to compute in the Lagrangian lattice formulation where one usu- 
ally computes the expectation value of the action. Finite density QED and QCD in the 
Lagrangian lattice formulation is hampered by the notorious complex action problem. 

- Atomic physics: study of spectra and the question of quantum chaos. 

- Condensed matter physics: study of spin systems (computation of dynamical structure 
factors), and high Tc superconductivity models (search for electron pair attraction at very 
small energy). In the following we will outline how to construct such a Hamiltonian. 



2 Construction of H^ff 

In contrast to the statistical mechanics concept of the transfer matrix, which describes 
the time-evolution (we consider imaginary time) when advancing the system by a small 
discrete time step At — qq and from which one can infer the Hamiltonian (ao — > 0), here 
we consider transition amplitudes < ■0|e~^^/'*|0 > corresponding to a finite, long time T 
{T » ao), for the purpose to reconstruct the spectrum in some finite low energy domain. 
Let us start from a complete orthonormal basis of Hilbert states |ej >, i = 1,2,3 - • • and 
consider the matrix elements for a given fixed N 

Mij{T) e,|e-^'^/'^|e,- >, j G 1, • ■ ■ , TV. (3) 

Under the assumption that H is Hermitian, M{T) is a positive, Hermitian matrix. Ele- 
mentary linear algebra implies that there is a unitary matrix U and a real, diagonal matrix 
D such that 

M(T) = D{T) U. (4) 

On the other hand, projecting H onto the the subspace Sn generated by the first N states 
of the basis |ei >, and using the eigenrepresentation of such Hamiltonian, one has 

M^AT) = E < e.\Ef > e-<^^/^ < Ef\e, >, (5) 
fe=i 
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and we can identify 

Ul=<e.\Elff>. D,{T) = e-K"^'\ (6) 

Let us assume for the moment that the matrix elements Mjj (T) , i, j = 1 , ■ ■ ■ , would be 
known. Then algebraic diagonalization of the matrix M(T) yields eigenvalues Dk{T), k = 
1, - ■ ■ ,N, which by Eq.(^) gives the spectrum of energies, 

El" = ~lnD,iT), k = l,...,N. (7) 

The corresponding k-th eigenvector can be identified with the k-th column of the matrix 
Uji.. From Eq.(P) we then know the wave function of the k-th eigenstate expressed in 
terms of the basis |ej >. Thus starting from the matrix elements Mij{T) we have explicitly 
constructed an effective Hamiltonian 

N 

Heff = j:m">Ei"<Ei"\. (8) 

k=l 

3 Computation of matrix elements by Monte Carlo 

We suggest to compute the matrix elements Mij{T) directly from the action via Monte 
Carlo with importance sampling. For the sake of simplicity, let us consider D = 1. We 
choose basis states |ej > in position space by introducing a lattice with nodes Xi and define 
ej(x) (unnormalized) by ej(x) = 1 if Xj < x < Xj+i, zero else. Axj = Xj+i — Xj. In numerical 
calculations we have used a regular lattice, Axi = const. The matrix elements read 



Mij{T) = I dyj dz<y,T\z,0> 



dy / dz [dx]exp[—S[x]/h] 



y,T 

(9) 

z,0 



Here 5* denotes the Euclidean action for a given path C, 



S\C] = r dt -mx^ + V{x) 
Jo 2 



(10) 

c 
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The Monte Carlo method with importance samphng is suited and conventionally applied 
to estimate a ratio of integrals, like in Eq.(|l|). Here we suggest to estimate the matrix 
elements M^- by splitting the action 



S — S(] -\- S\/ 



dt -mx^ + dt V{x), 
2 ^0 



(11) 



and to express Mij as 



M,,(T) = il4°)(T) 



Ix-^' dy J!^^+^ dz J[dx] exp[-Sv[x]/h] exp[-^oN/^] 



y,T 

z,0 



y,T 



z,0 



(12) 



where O = exp[—Sv/K\ is treated as an observable. The ratio can be treated by standard 
Monte Carlo methods with importance sampling. The matrix elements M-j^ , corresponding 
to the free action Sq, are almost known analytically, 



dy 



dz 



m 



27ihT 



exp 



m 



2hT 



{y - 



(13) 



4 Test oi He ff 
4.1 Free system 

In order to test the effective Hamiltonian, we have computed the energy spectrum, its wave 
functions and thermodynamic observables like the average energy U and the specific heat 
C as well as the partition function Z. They are defined by 



C{f3) 



dF 



logZ 



(14) 



dT '^"'^ 9/52 ' 

where jS = {kBT)~^, T is the temperature, and we identify (3 with the imaginary time T 
hj (3 = T/h. For the free system one obtains the following analytical expressions for Z, U 
and C, 



m 



dx (being infinite). 



CW) = (15) 

Note that U{P) — ^/j^oo 0, i.e. it tends to the ground state energy of the free system 
(Feynman-Kac formula). 

The partition function corresponding to the effective Hamiltonian is obtained via its 
spectrum, 

Z,ff{(3) = Trie-f'^^ff] = e"^^^''. (16) 

k=l 

Via Eq.(p^) one obtains the corresponding average energy Ueff and the specific heat Cg//. 
One should keep in mind that Hgff has been constructed for a specific value of the time 
parameter, T = 1 corresponding to the temperature T = 1 (we use h = ks = 1). Fig.[l] 
shows a plot of the average energy, comparing the exact result with that from the effective 
Hamiltonian. One observes that the agreement is better where T 0, i.e. in the low 
energy regime. A similar behavior is found for the specific heat, shown in Fig. [2]. 

4.2 Harmonic oscillator 

The Euclidean action of the harmonic oscillator is given by 

fT I 1 

S = dt—mx^ H — muj'^x^ . (17) 

Jo 2 2 

The energy spectrum is 

En = nuj{n + 1/2), n = 0,l,2,--- (18) 

A comparison of the spectrum of the effective Hamiltonian with the exact one is shown in 
Tab.[l]. As can be seen, the error is small in the low energy domain. A more stringent test 
is that of the wave functions. Fig. [3] shows a comparison for the wave functions of the three 
lowest states. We have also verified the low energy behavior of the effective Hamiltonian 
by computing the partition function, average energy and specific heat as a function of 



6 



temperature. For the harmonic oscillator those are analytically known, 



U{(3) 



[2 smh{phuj/2)]-^ 
^ ctgh{/3nLu/2) , 

r /3huj/2 



cm 



smh{phu/2) 



(19) 



In the limit /? — oo the average energy tends to the ground state energy, U huj/2 
(Feynman-Kac formula). A plot of the average energy and the specific heat is shown in 
Figs. [4,5]. The effective Hamiltonian, constructed at Tc = Pc = % = ^, describes well 
thermodynamic observables in the range I3c < {3 (it works also for (3 > 10, not shown in 
the figure). However, it breaks down for (3 < (3^, i-e. T > This is due to the small 
dimension = 20 of the matrix. Agreement in a larger /3-region, i.e. lowering (3c can 
be obtained by increasing A^. This can be seen, e.g. for the free system in Fig. [2], where 
N = 200 and < 0.1. 

4.3 Other local potentials 

We have tested the effective Hamiltonian for other local potentials. For example. 



is a potential having a minimum —Vq at x = and rising asymptotically to zero at x = ±oo. 
It generates a bound state spectrum, being analytically known 0. It is given by 



V{x) = -Vo sech.\x/d) 



(20) 



'n 



2md? 



X 



(n + 1/2) 



+ 1/4 



Q 



2m(P 



n 




(21) 



The results are shown in Tabs. [2,3]. 
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5 Conclusion 



We have proposed to construct an effective low-energy Hamiltonian from the action via 
use of the Monte Carlo method. We have shown that the method works for a number of 

systems in 1 — D quantum mechanics, by computing the spectrum, wave functions and 
thermodynamical observables. 

- We have not given an error estimate of the statistical errors. The reason is that the 
statistical error of the matrix elements can be estimated easily, however, to get from that 
an error estimate of the energy spectrum is difficult. We defer that to a later study. 

- We have not discussed an application to a field theory or a many-body system, although 
this is the area where the method should prove to be most useful. The reason is that 
this requires a new step, namely a stochastic (Monte Carlo) selection from the set of basis 
functions. This is presently under investigation. 

- In our opinion an effective low-energy Hamiltonian will be very useful in condensed matter 
physics, atomic physics, nuclear physics, and high energy particle physics. 
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Figure Caption 

Fig.l Average energy of the free system. Solid line and diamonds, respectively, represent 
the exact analytical result, and that from the exact matrix elements for Ax = 0.5 
and = 100. The cross at /5 = 0.1 corresponds to Ax = 0.2 and = 200. 

Fig. 2 Specific heat over ks of the free system. Symbols as in Fig. [1]. 

Fig. 3 Wave function of the harmonic oscillator, (a) ground state, (b) first excited state, 
(c) second excited state. Solid hne, diamonds and crosses, respectively, represent the 
exact analytical result, that from the exact matrix elements, and that from Monte 
Carlo simulation. 

Fig. 4 Average energy of the harmonic oscillator. Symbols as in Fig. [3]. 

Fig. 5 Specific heat over of the harmonic oscillator. Symbols as in Fig. [3]. 

Tab.l Eigenvalues of the harmonic oscillator, m = 1, h = 1, uj = 0.6, Ax = 1, N = 20. 
^exact^ ^e.m. ^^^^ -E'^''^', respectively, represent the exact analytical result, that from 
the exact matrix elements, and that from Monte Carlo simulation. 

Tab. 2 Bound state spectrum for potential given by Eq.(|20|). (a) For m = 1.0, h = 
1.0, T = 1.0, Vo = 1.0, d = 1.0, Q = 2, Ax = 1.0, A^ = 10, there is only 
one bound state Umax < 1- This is confirmed by the Monte Carlo data, (b) For 
m = 1.0, h = 1.0, T = 1.0, Vo = 1.0, d = 2.0, Q = 8, Ax = 1.0, A^ = 20, there are 
three bound states Umax < 3. This is confirmed by the Monte Carlo data. 
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Fig. [3b] 
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Fig. [4] 
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Harmonic oscillator 



m = l, h = l, uj = 0.6, Ax = 1, N = 20 
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Fig. [5] 
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